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Abstract 


This  report  an  effort  to  collect  under  one  cover  the  existing  theory  of 
three  continuous  zone  refining  systems;  the  zone-void  system,  the  zone-transport 
system  and  the  matter  transport  system.  The  steady- state  equations  which  de¬ 
scribe  these  systems  are  derived  using  consistent  terminology  and  under  boundary 
conditions  which  account  for  the  density  change  due  to  melting.  The  newly  devel¬ 
oped  nondime nsional  steady- state  equations  are  also  derived  and  the  results  of  an 
extensive  computer  study  based  onthese  equations  are  presented  in  graphical  form. 
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Continuous  Zone  Refining 


1.  INTRODLCTION 

That  trace  impurities  are  capable  of  influencing  the  physical  and  chemical 
properties  of  their  host  material  has  been  known  for  many  years.  Recently  in 
fact,  instances  have  been  found  in  which  trace  impurities  are  responsible  for  or 
control  critical  properties  of  the  matrix  material.  One  of  the  most  dramatic  ex¬ 
amples  of  the  ability  of  trace  impurities  to  control  a  property  of  its  matrix  is 
found  in  the  field  of  semiconductors.  The  variation  of  the  type,  quantity  or  dis¬ 
tribution  of  impurities  in  the  part  per  million  to  the  part  per  billion  range  can 
completely  cl  \nge  the  electriceil  properties  of  a  semiconductor.  For  this  reason 
the  control  of  trace  impurities  in  semiconductors  is  of  primary  importaix:e. 

In  order  to  systematically  study  the  effect  of  various  impurities  aiui  combina¬ 
tions  of  impurities  on  the  electrical  properties  of  a  semiconductor,  it  became 
necessary  to  prepare  semiconducting  materials  which  were  pure  to  the  part  per 
billion  range  or  better.  Tn  the  early  days  of  the  semiconductor  industry  the  most 
important  semiconductor  was  germanium.  Known  purification  technique.H  failed 
to  purify  germanium  to  the  required  purity  level.  Pure  germanium  was  finally 
prepared  by  a  new  purification  technique  developed  primarily  to  purify  semicon¬ 
ductor  materials.  This  technique,  which  was  originated  by  W.  G.  Pfann,  is  called 
zone  melting.  When  materials  are  purified  by  zone  melting  techniques  the  process 
is  termed  zone  refining. 

(Received  for  publication  22 October  1964) 


Zone  refining  is  not  a  technique  which  resulted  from  the  discovery  of  new 
principles.  In  retrospect  one  can  clearly  see  that  zone  refining  was  a  logical  ex¬ 
tension  of  the  use  of  the  freezing  or  crystallization  process  as  a  means  of  purifi¬ 
cation.  To  use  crystallization  as  a  means  of  purification  one  takes  advantage  of 
a  well  known  physical  fact,  namely,  when  a  solution  starts  to  freeze  the  solute  is 
almost  always  unevenly  distributed  between  the  liquid  and  solid  phases.  A  typical 
example  of  this  uneven  distribution  of  solute  is  found  in  the  use  of  the  freezing 
process  to  remove  salt  from  sea  water.  When  the  temperature  of  a  saline  solution 
is  lowered  below  the  freezing  point  the  first  ice  particles  which  freeze  out  have  a 
lower  concentration  of  salt  than  the  original  solution.  A  process  based  on  this 
fact  is  being  given  serious  consideration  as  a  means  of  the  large  scale  preparation 
of  fresh  water  from  sea  water  in  Israel. 

Under  equilibrium  condititms  the  measure  of  the  inequality  of  the  solute  dis¬ 
tribution  between  the  solid  and  liquid  phase  of  such  a  system  is  given  by  the  equi¬ 
librium  distribution  coefficient,  k^,  which  is  the  ratio  of  the  concentration  of 
solute  in  the  solid  Cg  to  the  solute  concentration  in  the  liquid  Cj^.  Thus 


(1) 


The  freezing  process  has  been  known  as  a  purification  technique  for  a  great 
many  years.  Its  utilization  however  has  been  very  limited  due  to  the  fact  that  it 
is  9X1  extremely  tedious  technique  to  employ.  Prior  to  zone  refining  the  only  basic 
improvements  made  in  art  of  purification  by  freezing  was  fractional  crystalliza¬ 
tion  and  fractional  crystallization  was  simply  a  repetition  of  the  freezing  process. 

In  fractional  crystallization  just  as  in  a  regular  crystallization  all  of  a  solution 
consisting  of  two  or  more  components  is  made  molten.  The  solution  is  then  allowed 
to  cool  until  a  predetermined  fraction  of  the  total  solution  has  solidified.  The  solid- 
fied  fraction  is  then  separated  from  the  melt  by  mechanical  means.  In  a  simple 
crystallization  (assuming  a  favorable  k^  and  an  ideal  binr^ry  mixture)  the  partial 
separation  of  the  two  components  obtained  in  this  manner  constitutes  the  entire 
purification  procedure.  In  a  fractional  crystallization  this  process  would  be  re¬ 
peated  using  the  separated  solidified  fraction  of  the  first  crystallization  as  the 
starting  solution  of  the  second  crystallization.  The  number  of  crystallizations  used 
would  depend  on  the  purity  of  the  product  desired.  Probably  the  most  famous  ap¬ 
plication  of  fractional  crystallization  was  the  separation  of  radium  from  barium  by 
the  Curies.  The  fact  that  by  use  of  this  technique  the  Curies  and  others  were  able 
to  effect  the  remarkable  separations  they  did  is  a  tribute  to  their  analytical  tech¬ 
nique  and  their  determination  rather  than  the  efficiency  of  the  technique. 
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Purity  via  fractional  crystallization  is  an  arduous,  tedious  chore.  This  process 
requires  that  the  material  being  purified  imdergo  a  change  in  state.  The  drawback 
of  the  process  is  that  one  of  these  states  is  the  solid  state.  The  reason  tliat  the 
solid  state  poses  such  difficulties  in  a  purification  process  is  that  under  normal 
conditions  solids  do  not  flow.  When  the  solid  state  is  met  in  such  a  process,  the 
separation  of  pure  and  impure  fractions  must  be  effected  by  the  operator  and  not 
by  the  materials  themselves  (as  is  the  case  in  distillation).  Unfortunately,  until 
zone  refining  was  discovered  fractional  crystallization  was  the  only  purification 
method  employing  the  freezing  principle,  by  which  a  separation  factor  of  any  given 
multiple  of  the  distribution  coefficient  could  be  obtained. 

The  production  of  single  crystals  from  the  melt  is  an  application  of  the  freez¬ 
ing  process.  Therefore  under  ideal  conditions  for  crystal  growth,  the  measure 
of  the  distribution  of  impurities  between  the  solid  and  the  liquid  at  the  freezing 
interface  is  given  by  the  distribution  coefficient.  In  the  Bridgeman^  process  for 
crystal  growth,  where  both  the  melt  and  the  growing  crystal  are  contained  in  the 
same  crucible  and  solidification  proceeds  slowly  from  one  end;  if  the  distribution 
coefficient,  k,  is  less  than  one,  the  growing  crystal  injects  impurities  into  the 
liquid.  When  crystallization  is  complete  the  impurities  will  be  unevenly  distributed 
throughout  the  crystal. 

Making  the  following  assumptions: 

1.  The  original  solution  is  an  ideal  binary  mixture 

2.  Diffusion  in  the  solid  is  negligible 

3.  Uniform  solute  distribution  in  the  liquid 

4.  k  (the  effective  distribution  coefficient  of  the  growing  crystal)  is  constant 

2 

Pfann  r^orts  that  C,  the  concentration  of  impurity  at  any  plane  along  the  crystal, 
can  be  expressed  by  Eq.  (2).  Thus, 

C  =  kC^d-g)'^"^  (2) 

where 

C  =  solute  concentration  in  the  solid  at  the  liquid-solid  interface,  in 
units  of  solute  per  unit  volume  of  solu';ion. 

=  initial  impurity  concentration 

and 

g  =  the  fraction  of  the  original  solution  that  has  solidified. 

2 

Pfann  called  this  type  of  freezing  process  normal  freezing.  Equation  (2)  of  course 
assumes  a  crystal  or  ingot  of  uniform  cross  section  whereas  in  the  Bridgeman 
technique  the  ingot  usually  tapers  to  a  point  at  one  end. 


wimif 
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In  19u?;  Pfann  reported  the  basic  principles  of  zone  melting.  Zone  melting 
is  a  technique  whereby  the  freezing  process  could  be  carried  out  on  an  ingot  in 
such  a  way  that  impurities  could  be  distributed  in  a  manner  which  differed  from 
that  which  would  be  produced  by  a  normal  freeze.  The  main  difference  between 
zone  melting,  fractional  crystallization  and  normal  freezing  is  in  the  amount  of 
material  that  is  made  molten  at  any  one  time.  In  both  fractional  crystallization 
and  normal  freezing  all  of  the  material  undergoing  treatment  is  made  molten  at 
the  start  of  the  process.  For  zone  melting  only  a  fraction  of  the  total  material  is 
made  molten  at  any  one  time.  Zone  melting  is  nvjrmally  accomplished,  as  shown 
in  Figure  1,  by  arranging  the  material  to  be  treated  in  the  shape  of  a  bar  and 
causing  a  molten  zone  to  travel  from  one  end  to  the  other.  The  solute  distribution 
which  results  from  the  passage  of  such  molten  zones  along  a  bar  is  a  function  of 
many  parameters  among  which  are:  the  distribution  coefficient,  the  original  solute 
distribution,  the  number  of  zones  which  traverse  the  bar,  the  direction  in  which 
these  zones  travel,  the  size  of  the  zones  and  the  length  of  the  bar.  It  is  obvious 
that  a  great  many  different  solute  distributions  are  possible.  The  methods  for  con¬ 
trolling  these  solute  distributions  by  variation  of  these  parameters  are  called  zone 
melting  techniques. 

The  most  important  application  of  zone  melting  is  in  the  purification  of  mate¬ 
rials.  Zone  melting  techniques  which  are  designed  specifically  to  purify  materials 
are  called  zone  refining  techniques.  For  simplicity  the  discussions  of  zone  re¬ 
fining  which  follow  will  be  restricted  to  the  purification  of  an  ideal  binary  solution 
in  which  the  solute  (impurity)  is  present  in  very  low  concentrations  and  k  is  con¬ 
stant.  In  this  system  we  also  assume  that  there  is  no  diffusion  in  the  solid  and 
complete  diffusion  in  the  liquid.  These  assumptions  which  will  hold  for  the  re¬ 
mainder  of  this  report  are  very  closely  approximated  in  most  zone  refining  systems. 

The  equilibrium  distribution  coefficient  can  be  determined  from  the  phase 
diagram  of  the  binary  system  being  purified.  If  the  solute  raises  the  melting  point 
of  the  solvent  then  k^  >  1;  however,  if  the  solute  lowers  the  melting  point  of  the 
solvent  then  k^  <  1.  The  value  of  k^  can  be  attained  only  under  perfect  equilibrium 
conditions.  In  an  actual  zone  refining  operation  the  value  of  k^  can  only  be  ap¬ 
proached.  The  value  of  the  effective  distribution  coefficient,  k,  the  distribution 
coefficient  which  actually  exists,  can  and  must  be  determined  experimentally  for 
any  given  system.  If  normal  experimental  care  is  exercised  the  assumption  of  a 
constant  k  for  dilute  solutions  is  quite  reasonable. 

Purification  via  zone  refining  is  accomplished  by  passing  a  number  of  molten 
zones  through  an  ingot  of  the  material  to  be  purified.  All  the  zones  travel  in  the 
same  direction.  The  number  of  zones  is  predetermined  by  the  purity  desired. 

The  direction  of  travel  of  impurities  is  a  function  of  the  effective  distribution  co¬ 
efficient.  For  k  <  1,  impurities  collect  in  the  liquid  and  travel  with  the  zones. 
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For  k  >  1,  impurities  are  rejected  into  the  solid  and  travel  in  the  direction  opposite 
to  the  zones  travel.  The  degree  of  purification  is  a  function  sf  the  number  of  zones 
passed  through  the  ingot  and  approaches  a  limit  as  the  niunber  of  zones  used  ap¬ 
proaches  infinity. 

The  equation  which  describes  the  solute  distribution  in  an  ingot  after  one  mol- 

2 

ten  zone  (pass)  has  been  put  through  the  ingot  was  reported  by  Pfann.  This  equa¬ 
tion  Pfann  expressed  as: 

kx 

^  =  1  -  (1-k)  e  (3) 


Equation  (3)  describes  the  solute  distribution  which  exists  in  all  but  the  last  zone 
length  of  the  ingot.  A  complete  derivation  of  Eq.  (3),  the  idealizing  assumptions 
on  which  it  is  based,  and  the  definitions  of  the  terms  involved  are  given  in  Appen¬ 
dix  A.  An  equation  describing  the  solute  distribution  in  an  ingot  following  the  pas- 

3  4 

sage  of  more  than  one  molten  zone  is  not  easily  derived.  Lord  and  Reiss  were 

5 

the  first  to  present  approximate  solutions.  It  was  left  to  Braun  and  Marshall  to 
generate  a  complete  analytical  solution  for  multipass  solute  distributions.  Even 
with  Braun's  expressions,  however,  obtaining  the  solute  distribution  which  exists 
after  five  or  more  passes  requires  a  great  deal  of  computation  and  is  a  problem 
which  is  more  realistically  solved  by  a  computer. 

The  separation  of  impurities  obtainable  by  passing  molten  zones  through  a 
charge  increases  as  the  number  of  passes  increases  until  a  steady- state  or  ulti¬ 
mate  distribution  is  reached.  When  the  steady-state  is  obtained,  the  passage  of 

additional  zones  through  the  ingot  does  not  alter  the  solute  distribution.  The  equa- 

2 

tions  which  describe  this  distribution  were  first  presented  by  PfEuin.  They  are: 


C(x)  =  Ae®* 
where 


k  = 


B£ 

B£ 

e 


and 


A  = 


C  BL 
o 

BL  , 
e  -  1 


(4) 


(5) 


(6) 
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A  derivation  of  these  equations  is  given  in  Appendix  B.  These  equations  do  nol 
take  into  account  the  back  reflection  of  solute  due  to  the  effect  of  the  normal  freeze 

g 

in  the  terminal  zone  of  the  charge.  Braun  presents  a  solution  for  the  ultimate 
distribution,  for  any  value  of  k,  which  takes  into  account  the  effect  on  the  ultimate 
distribution  due  to  the  normal  freeze  which  occurs  in  the  final  zone  during  each 
pass. 

Zone  refining  is  an  easy  way  to  carry  out  a  fractional  crystallization.  The 
simple  change  in  procedure  namely,  the  passage  of  molten  zones  has  greatly  in¬ 
creased  the  efficiency  of  crystallization  as  a  purification  technique  in  that  it  has 
automated  the  separation  of  pure  and  impure  fractions.  A  number  of  molten  zones 
can  be  passed  through  an  ingot  in  a  single  operation  providing,  in  simple  fashion, 
the  end  of  fractional  crystallization  which,  as  has  already  been  mentioned,  is  a 
separation  factor  which  is  any  given  multiple  of  the  distribution  coefficient. 


2.  CONTINIJOLS  ZONE  REFINING 

Zone  refining  as  described  by  Pfann  in  Reference  2  is  a  multistage  batch  pro¬ 
cess.  Such  a  batch  process  can  in  many  instances  fulfill  all  of  the  requirements 
involved  in  the  production  of  a  pure  material.  Batch  processes,  however,  have 
many  intrinsic  limitations  which  can  be  eliminated  if  the  process  can  be  made  con¬ 
tinuous.  In  general,  a  continuous  purification  process  is  one  which  comprises  an 
enriching  and  a  stripping  section  and  the  means  for  the  introduction  of  feed  and 
the  removal  of  product  and  waste.  The  addition  of  feed  to  a  continuous  system 
must  occur  at  a  point  in  the  system  which  is  between  the  places  where  product  and 
waste  are  removed.  Continuous  fractional  distillation  is  the  analog  of  continuous 
zone  refining  which  treats  the  phase  change  between  vapor  and  liquid.  In  principle 
they  are  very  similar.  In  practice  the  ability  of  the  vapor  phase  to  be  compressed 
easily  and  the  fact  that  the  liquid  phase  flows  and  does  not  build  up  on  the  contziiner 
walls  allows  the  feed  to  be  introduced  easily  at  any  point  on  a  distillation  column 
and  also  provides  counter  current  flow.  In  a  zone  refining  system  things  are  not 
quite  so  simple.  Here  neither  phase  is  easily  compressed  and  the  solid  instead  of 
flowing  freely  usually  adheres  tenaciously  to  and  builds  up  on  the  container  walls. 
The  continuous  zone  refiner  must  live  with  these  difficulties  and  still  provide  for 
the  requisite  flows  of  feed,  product  and  waste  while  permitting  the  unidirectional 
movement  of  molten  zones  down  the  column.  Such  a  system  is  shown  in  generalized 
form  in  Figure  2. 

The  creation  of  molten  zones  and  the  mechanism  by  which  these  zones  are 
made  to  traverse  the  column  are  experimental  problems  and  therefore  are  beyond 
the  scope  of  this  paper.  It  will  be  assumed  that  the  mechanism  for  their  creation 


Figure  1.  A  Molten  Zone  Being  Passed  Through  an  Ingot  of  Length  L 


DIRECTION  OF  TRAVEL  OF  HEATERS 


Figure  2.  A  General  Continuous  Zone  Refining  Apparatus 
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and  movement  has  been  found  and  that  this  mechanism  produces  zones  which  have 
planar  interfaces,  are  of  constant  size  and  move  at  a  rate  of  travel  which  is  con¬ 
stant. 

7 

In  1955,  Pfann  described  the  zone-void  system,  the  first  zone  refining  system 
postulated  capable  of  completely  continuous  operation.  In  this  system  the  movement 
of  material  into  and  out  of  a  system  is  brought  about  by  the  movement  of  voids  in¬ 
side  of  the  system.  Voids  are  normally  gas  filled  spaces,  however,  in  a  more 

g 

general  fashion  Pfann  has  defined  a  void  as  a  region  in  which  there  is  a  substantial 

deficiency  of  the  material  undergoing  treatment.  Material  movement  is  produced 

by  causing  voids  formed  at  the  product  and  waste  exits  to  move  to  the  feed  inlet. 

The  column  is  constructed  in  two  parallel  vertical  sections  joined  together  at  the 

top  by  the  feed  section.  An  example  of  this  type  of  apparatus  is  shown  in  Figure  3. 

Voids  are  created  at  the  bottom  of  each  section  by  the  controlled  removal  of  a 

small  amount  of  material  from  the  colamn.  The  size  of  the  void  must  be  less 

8  9 

than  one  molten  zone.  Several  methods  of  forming  voids  have  been  reported.  ' 
Product  and  waste  are  collected  at  the  bottom  of  the  enriching  and  stripping  sec¬ 
tions,  respectively,  in  forming  the  voids.  Material  is  moved  down  the  column  as 
a  result  of  movement  of  the  molten  zones  which  cause  the  voids  to  move  upward. 
When  the  voids  reach  the  top  of  the  section  they  are  filled  in  by  feed.  Thus,  feed 
is  introduced  into  a  column  every  time  a  void  reaches  the  feed  section.  In  the 
enriching  section,  voids  travel  continuously  with  the  molten  zones  while  in  the 
stripping  section  their  travel  is  intermittent,  a  zone  length  at  a  time. 

One  of  the  primary  advantages  of  continuous  operation  is  the  ability  to  operate 
at  or  near  the  steady  state.  The  time  required  to  attain  this  distribution  will  vary 
from  system  to  system,  however,  even  if  the  start  up  time  is  long,  once  atteiined, 
the  system  continues  to  deliver  material  of  maximum  purity. 
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Pfann  has  shown  that  for  any  continuous  zone  refining  system  the  overall 
material  balance  in  terms  of  the  volume  flow  rates  of  feed,  F,  waste,  W,  auid 
product,  P,  is  given  by: 

F  =  W  +  P  (7) 

Ii\  a  continuous  zone  refining  system  it  is  the  function  of  tlie  enriching  section  to 
produce  a  pure  product.  In  doing  this  the  enriching  section  injects  into  the  feed 
the  impurities  which  it  has  removed  from  the  product  it  produces.  Without  a 
stripping  section  this  injection  of  impurities  into  the  feed  would  eventually  increase 
the  solute  concentration  in  the  feed  to  the  point  where  no  purification  wo  ,’.d  be  ob¬ 
tained.  The  function  of  the  stripping  section  is  to  inject  into  the  feed  material  of 
sufficient  purity  and  in  sufficient  quantity  to  maintain  the  solute  concentration  in 
the  feed  constant.  When  the  stripping  section  accomplishes  this  function  the  solute 
balance  for  the  system  is  given  by 


10 


FCj  =  PCp  +  WC^  (8) 

where  C-,  C  ,  and  C  are  the  solute  concentrations  in,  respectively,  the  feed, 

1  p  w 

product,  and  waste.  When 

Cp  =  aCj  (9) 

and 

=  PCj  (10) 

the  ratio  of  the  flows  of  product  to  waste  (P/W),  which  will  provide  this  solute 
balance  can  be  found  by  substituting  Eqs.  (7),  (9)  and  (10)  into  Eq.  (8).  Thus 


P. 

W 


1  -  a 


(11) 


Purity  requirements  and  the  solute  concentration  in  the  feed  determine  a. 
When  the  steady  state  has  been  reached  the  weight  of  solid  material,  L^,  that  must 
be  in  the  enriching  section  to  effect  the  constant  removal  of  a  product  of  concen¬ 
tration  Cp  can  be  obtained  by  solving  the  steady- state  equations.  Thus,  for  a 
material  of  known  k,  in  a  system  using  a  molten  zone  of  convenient  weight,  I,  and 
a  convenient  void  size,  is  given  by 

^e  0.43^  Bg  [^^(4  ■  O]  <^2) 

where  B^,  w  and  i//  are  constants  which  depend  on  h,  jf  and  k  and  are  given  by: 


,(h-l)' 


(14) 


.b  -  hk-£ 
-  TTT 


(15) 


n 

The  derivation  of  Eqs.  (13)  to  (15)  was  first  reported  by  Pfann.  The  authors  have 
modified  the  definitions  of  terms  originally  used  by  Pfuin  to  account  for  the  change 


11 


in  density  due  to  melting.  The  derivation  of  Eqs.  (13)  to  (15)  along  with  the  def¬ 
initions  of  previously  undefined  terms  is  given  in  Appendix  C. 

The  steady-state  solute  distribution  these  equations  describe  is  shown  in 
Figure  4.  The  mrithematical  model  does  not  give  any  information  about  the  actual 

7 

total  weight  of  material  in  the  enriching  section,  however,  Pfann  re[}orts  that  in 

practice  this  weight  is  usually  equal  to  L  +  h. 

7  ® 

Pfann  also  presented  the  steady-state  equations  for  the  stripping  section  of 

the  zone- void  system.  As  a  result  of  his  derivation  Pfann  reported  that  in  the 

stripping  section  the  length  of  column  L  which  will,  produce  a  separation  P  is 

s 

given  by 


S  '  oiwb;  ‘“s  [(^)  *  ‘] 

where  B  ,  6  and  y  are  constants  and  functions  of  h,  £ ,  k  and  p  which  can  be 
s 

obtained  by  the  following; 

‘ 

[(h-£)p-h(l-k)] 

y  «  (19) 

1  +  B  h-e  ® 
s 

As  part  of  the  theoretical  continuation  of  the  zone  refining  effort  in  this  labo¬ 
ratory  a  study  of  the  effect  of  variation  of  the  significant  parameters  of  the  zone- 
void  system  was  undertaken.  The  parameters  under  study  were  h,  £,  and  k  (as 

7 

defined  by  Pfann  ).  The  object  of  this  study  was  to  determine  the  effect  of  these 
parameters  on  a  and  p. 

in  the  stripping  section  the  definition  of  k  and  p  precludes  the  possibility 
of  p  assuming  values  less  than  unity  when  k  <  1.  In  the  course  of  this  investi¬ 
gation  it  was  found  that  computations  based  on  Eqs.  (16)  to  (19)  gave  values  of  p 

less  th£ui  unity  in  the  region  1  >  k  >  £/h  for  almost  all  values  of  L  .  Further 

s 

computations  showed  that  p  assumes  negative  values  when  k  >  1  for  many  values 

of  L  ,  £  and  h.  A  study  of  Pfann's  derivation  of  Eqs.  (16)  to  (19)  revealed  an 
s 

error  in  the  equation  describing  the  solute  balance  over  the  entire  cycle. 


(17) 

(18) 


I 

i 


f 

j 

I 


I 

I 
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(Equation  (28),  Reference  7.)  This  error  invalidated  tlie  final  equations  of  Pfann's 
derivation. 

The  equation  cf  the  correct  solute  balance  is 

J 

hCj=  /  C(x)dx+ (h-l)p  Cj  (20) 

o 

Using  Eq.  (20),  corrected  steady-state  equations  which  account  for  the  change 
in  density  due  to  melting  are  derived  as  shown  in  Appendix  D.  In  this  derivation 
terms  are  redefined  to  account  for  the  density  change.  Based  on  these  corrections 
and  the  accompanying  definition  changes  is  given  by 

K  =  OjW  ‘“8[fTT^]  '2‘' 


where 


B  (£-w)  , 

e  ®  1  f  (22) 


(23) 


and 


-  k(£-w)  (2^> 

Calculations  based  on  Eqs.  (21)  to  (24)  remove  the  aforementioned  inconsistencies 

in  p.  Figure  5  shows  the  steady-state  solute  distribution  in  the  stripping  section. 

A  paper  describing  part  of  this  phase  of  the  effort  has  been  published  elsewhere.^® 

After  the  error  in  Eq.  (28)  of  Reference  7  and  its  implications  with  respect 

11  12 

to  the  zone -void  system  had  been  noted  (as  above),  Pfann  pointed  out  that  Abe 

using  an  independent  method  had  derived  steady-state  equations  for  the  zone- 

transport  system  and  had  already  reported  the  existance  of  this  error  in  a  paper 

which  treated  the  zone-transport  system.  The  zone-transport  system,  a  continuous 

13  14 

zone  refining  system,  was  also  invented  by  Pfann.  '  An  example  of  a  zone- 
transport  apparatus  is  shown  in  Figure  6.  In  this  system  movement  of  material 
from  the  feed  section  to  the  product  and  waste  exits  is  accomplished  by  virtue  of 
the  difference  in  level  between  the  feed  and  the  product  and  waste  exits.  Because 


Sotute  Concentrotton.C,  in  Weight  Frociton 
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Figure  4.  Steady-State  Solute  Distribution  in  the  Enriching 
Section  of  the  Zone -Void  System 


Figure  5.  Steady-State  Solute  Distribution  in  the  Stripping  Section  of  the 
Zone-Void  System  when  the  Terminal  Zone  Undergoes  a  Normal  Freeze 


Wotle  Exit 
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of  this  difference  in  level,  material  in  the  molten  zones  flows  down  hill.  The 

weights  of  material  produced,  per  pass,  at  the  product  and  waste  exits  respectively 
14 

were  given  by  Pfann  as 


P  = 


v: 

2  L 


w  = 


2  L. 


(25) 


(26) 


In  describing  conditions  existing  when  the  steady  state  has  been  obtained, 

14 

Pfann  stated  that  the  steady-state  equations  which  describe  the  zone-void  system 

could  be  used  to  describe  the  zone-transport  system  provided  the  same  assumptions 

are  made  and  that  the  expressions  for  the  parameters  which  define  material  flow 

were  changed  in  accord  with  Table  1 .  CXir  investigation  has  shown  that  this  is  not 

a  valid  transition  due  to  the  fact  that  the  solute  concentrations  in  the  product  and 

waste  which  are  used  as  parameters  in  the  derivation  of  the  steady- state  equations 

are  different  in  the  zone-void  and  in  the  zone-transport  systems.  Thus  in  the  zone- 

void  system  C  is  given  by  Eq.  (27)  and  C  by  Eq.  (28). 
p  w 

TABLE  1 .  Expressions  for  the  material  flow  parameters 
in  each  system  when  constant  density  is  assumed 


Pfann 

Pfann 

Abe 

Zone- 

Zone- 

Zone- 

Void 

Transport 

Transport 

h  -  i 

P 

P 

Enriching 

t. 

£ 

£-  P 

Section 

h 

P  +  £ 

£ 

h  -  i 

W 

W 

Stripping 

h 

£ 

£  +  W 

Section 

£ 

£  -W 

£ 

Thus 


1  r''"* 


(27) 
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The  value  of  C^,  depends  upon  the  mode  of  freezing  in  the  last  zone.  In  the 

zone- void  system  is  defined  as  the  average  solute  concentration  of  the  charge 

between  L  +  £  and  L  +  h. 
s  s 

Thus 


w 


iT  I  f  ®  C  (x)  dx 
s 


(28) 


In  the  zone -transport  system,  because  of  the  physical  difference  in  the  mode  of 
collection  of  the  product  and  the  waste.  Cp  is  given  by  Eq.  (29). 


=  T  /  C(x)dx 


(29) 


while 


C 


w 


k 


(30) 
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In  his  study  of  the  zone-transport  system  Abe  transformed  Pfann's  zone- 

void  steady-state  equations  into  Pfann's  zone -transport  equations  by  following 

14 

Pfann's  directions.  Using  Pfann's  zone- transport  steady-stale  equations  Abe 

made  calculations  which  revealed  the  same  inconsistencies  in  P  noted  previously. 

12 

Using  his  own  technique  Abe  derived  steady-state  equations  for  the  zone- 

transport  system.  These  derivations  are  given  in  Appendices  E  and  F.  Using  his 
12 

own  equations  Abe  found  that  all  of  the  inconsistencies  in  q  and  p  were  elimi¬ 
nated. 

12 

In  his  paper  Abe  compared  the  limits  of  a  and  P  obtained  when  his  zone- 
transport  equations  and  Pfann's  zone-transport  '  quations  were  used  (at  constant 
density).  This  comparison  revealed  three  main  discrepancies.  These  discrep¬ 
ancies  are: 

1.  In  the  enriching  section  when  kf  1  and  ==  0,  Pfann  reports  q=  un/z/l  +0) 

whereas  Abe  finds  q  =  1  when  L  =  0  irrespective  of  the  value  of  k.  Where  in  the 

®  14 

zone  transport  system,  o),  L^.  £  and  P  are  defined  by  Pfann  as  follows: 
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k  + 


(k~l)£ 

P 


L  =  the  distance  from  the  exit  end  of  the  enriching  section  to  the  plane 
at  which  the  solute  cone,  in  the  solid  is  kC^ 

I  =  the  length  of  the  molten  zone 

P  =  the  volume  of  product  collected  per  pass 


2.  In  the  stripping  section  when  the  actual  length  of  the  stripping  section  minus 
one  zone  length  (Pfann  calls  this  length  L  .  Abe  calls  it  L*  )  equals  infinity  and 

k  =  1,  Pfann' s  equations  show  p  =  0  whereas  Abe's  equations  yield  p  =  1. 

3.  In  the  stripping  section  when  the  distance  defined  as  the  actual  length  of 

the  stripping  section  minus  one  zcxie  length  is  large  enough  and  when  k  >  1,  Pfann' s 

equations  yield  p  <  0  whereas  under  the  same  conditions,  Abe's  equations  have 

12 

p  =  0  when  k  >  (£+w)/f  (where  1  and  w  are  defined  by  Abe). 

15 

The  authors  have  pointed  out  that  in  their  opinion  these  differences  and  other 
differences,  which  become  apparent  when  References  7  and  12- 14  aire  studied  care¬ 
fully,  are  the  result  of  three  factors: 

1.  Partially  inconsistent  definition  of  terms. 

2.  The  error  in  Eq.  (28)  of  Reference  7. 

3.  The  transformation  of  the  steady-state  equations  of  the  zone-void  system 
to  the  steady-state  equations  of  the  zone-transport  system  without  taking  into  ac¬ 
count  the  difference  in  the  mathematical  expressions  of  the  product  and  waste  con¬ 
centrations  in  each  system. 

The  derivation  of  the  steady-state  equations  for  the  zone-transport  system  as 
14  12 

reported  by  Pfaxui  and  Abe  are  different  and  are  based  on  different  assumptions. 
However,  in  the  enriching  section  of  each  system  a  comparison  of  the  expressions 
that  define  materird  flows  is  straightforward  since  the  physical  and  mathematical 
significance  of  these  flows  is  identical.  For  the  stripping  section  a  comparison  of 
the  expressions  that  define  material  flows  is  complicated  by  the  different  methods 
of  derivation.  Thus,  in  the  stripping  section,  the  physical  volume  of  the  zone  that 
is  maintained  molten  and  caused  to  travel  by  a  moving  heater  is,  for  the  purpose 
of  each  derivation,  different.  The  mathematical  significance  of  the  molten  zone  in 
each  derivation,  however,  is  identical.  Table  1  lists  the  expressions  for  the  param¬ 
eters  that  define  the  saune  material  flows  in  Pfann's  zone-void  system  as  well  as 

Pfann' s  and  Abe's  zone-transport  systems.  The  differences  in  definition  of  terms 

12 

between  systems  which  Table  1  clearly  shows  are  not  explicit  in  Abe's  paper. 

The  definitions  of  terms  used  to  define  column  lengths  are  also  partially  in¬ 
consistent  or  different.  For  the  zone-void  system  (and  hence  the  zone-transport 
14 

system)  Pfann  defines  as  the  distance  from  the  exit  end  of  tlie  enriching 
section  to  the  plane  at  which  the  solute  concentration  in  the  solid  is  kC^  and  com¬ 
ments  that  in  practice  this  length  will  usually  be  one  heater  length  less  than  the 
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14 

actual  section  length.  For  the  zone-transport  system,  Pfann  defines  the  actual 
length  of  the  enriching  section  as  Lg.  In  Reference  12,  Appendix  B,  Abe  defines 
the  same  way  Pfann  does;  however,  in  Appendix  A  and  Figures  la  and  2  of 

Reference  12,  L  is  shown  as  the  actual  section  length. 

®  14 

In  the  stripping  section  of  the  zone-void  system,  Pfann  defines  the  length 

Lg  by 

S  H 

where  L*  is  the  actual  length  of  the  stripping  section  and  n  is  the  number  of  void 
encounters.  Since  there  are  no  voids  in  the  zone-trauisport  system,  Pfann' s  def¬ 
inition  of  L  for  the  zone-transport  system  is  given  by 

L  =  L'  -  h 
s  s 

Abe  on  the  other  hand  defines  L  as  the  actual  which  is  the  reverse  of  Pfann' s 

s 

terminology.  Thus  for  Abe 

K  -  h 
s  s 

The  authors  can  find  nothing  wrong  with  the  definition  of  terms  used  by  I  .ann  to 
describe  his  zone- void  or  zone-transport  systems.  Therefore  since  Pfann  invented 
both  systems  and  because  of  the  almost  universal  acceptance  of  his  terminology 
except  where  difficulties  were  encountered  in  treating  the  density  change  due  to 
melting,  the  authors  have  tried  to  adhere  rigidly  to  his  terms  and  definitions  of 
terms. 

This  is  true  even  in  Appendices  E  and  F  where  we  use  Abe's  derivation  of  the 
steady- state  equations. 

In  his  derivation  of  the  steady-state  equations  for  the  zone -transport  system 
Abe  introduces  the  term  y,  Pg/P£.  fo  account  for  the  density  change  due  to  melt¬ 
ing.  The  authors  feel  that  a  simpler  method  of  accounting  for  this  density  change 

in  all  continuous  zone- refining  systems  is  to  redefine  terms  by  expressing  solute 

1 6 

concentration  in  weight-fraction  as  suggested  by  Pfann  and  to  change  the  units 
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of  measure  of  position  along  a  column  from  length  to  weight  of  solid  material. 

When  these  changes  are  made  the  expressions  for  the  material  flow  parameters 
and  their  respective  definitions  for  each  system  are  given  in  Table  2. 

Abe's  derivation  of  the  steady-state  equations  for  the  zone-transport  system 
as  given  in  Appendices  E  and  F  have  been  modified  in  light  of  these  decisions  con¬ 
cerning  the  definition  and  consistency  of  terms  auid  the  handling  of  the  density 


TABLE  2.  Expreseions  for  the  material  flow  parameters 
in  each  system  when  the  density  change  due  to  melting  is 
taken  into  account 


Enriching 

Section 


Stripping 

Section 


Pfann 

Pfann 

Abe 

Zone- 

Zone- 

Zone- 

Void 

Transport 

Transport 

h  -  £ 

P 

P 

P  +  f 


1  -  P 


£  +  w 


Definition 


Weight  of  product  collected 
per  pass 

Weight  of  liquid  material 
in  the  molten  zone  that  is 
maintained  molten  and 
caused  to  travel  by  a 
moving  heater 

Weight  of  material  equal  to 
the  weight  of  liquid  mate¬ 
rial  that  is  maintained 
molten  and  caused  to  trav¬ 
el  by  a  moving  heater  plus 
the  weight  of  product  col¬ 
lected  per  pass 

Weight  of  waste  collected 
per  pass 

Weight  of  liquid  material 
in  the  molten  zone  that  is 
maintained  molten  and 
caused  to  travel  by  a  mov¬ 
ing  heater 

Weight  of  solid  material 
that  would  fill  the  volume 
occupied  by  the  molt <^n  zone 
when  the  cross  section  of 
the  column  is  unity.* 


*In  Abe's  zone-transprort  system  this  volume  is  smaller  than  the 
volume  of  the  molten  zone. 


change  due  to  melting.  These  changes  correct  an  error  in  the  handling  of  density 

present  in  Abe's  derivation  of  the  steady-state  equations  in  the  stripping  section.^^ 

A  derivation  of  the  steady-state  equations  for  the  enriching  section  of  the  zone- 

transport  system  is  given  in  Appendix  G.  This  is  the  same  derivation  outlined  by 
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Pfann  except  that  it  has  been  modified  in  light  of  Table  2  and  to  account  for  the 
different  product  concentrations  which  exist  in  the  zone-void  and  the  zone-transport 
systems.  As  a  result  of  this  modification,  the  expression  for  is  given  by; 


0.434  B 


log  .  u,l 

^  L  “(p+i)  ‘"J 
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where 


w  = 


P 

(P+l)  (1-k) 


and 


Bg(P+£) 


1  + 


(32) 


(33) 


In  the  stripping  sections  of  both  the  zone-void  and  zone-transport  systems  the 
shift  in  the  solute  distribution  due  to  the  mode  of  freezing  of  the  final  zone  at  the 
end  of  each  pass  has  not  been  treated  analytically.  Therefore  due  to  the  fact  that 
the  waste  is  withdrawn  from  the  zone-transport  system  before  the  last  zone  begins 
to  freeze  the  solute  concentration  in  the  waste  (back  reflection  not  considered)  is 
independent  of  the  mode  of  freezing  in  the  terminal  zone.  The  steady-state  equa¬ 
tions  for  the  stripping  section  of  the  zone-transport  s>  stem  can  be  derived  as  shown 
in  Appendix  H.  This  is  the  same  derivation  used  by  Pfann  to  describe  the  stripping 
section  of  the  zone-void  system  corrected  for  the  error  in  Eq.  (28)  of  Reference  Y 

and  using  the  proper  waste  concentration  (Eq.  (30)).  Using  this  derivation  L  is 

s 

given  by 

^  ‘-e  [^] 


where 


B  (£-w)  B  £ 

e  ®  =  1  +  -f- 

k 


(35) 


and 


'^"k(w-jr)+£ 

Abe's  derivation  of  the  zone-transport  equations  (Appendices  E  and  F)  differs 
substantially  from  Pfann' s  corrected  derivation  (Appendices  G  and  H).  However, 
both  derivations  describe  the  same  system  and  unless  the  mathematical  model  is 
different  the  expressions  for  which  result  should  be  identical.  Algebraic  ma¬ 
nipulation  using  Table  2  shows  that  this  is,  indeed,  the  case. 


The  iimits  of  a  and  p  therefore  are  also  identical.  These  limits  and  those 
for  the  zone  void  system  are  shown  in  Table  3. 

TABLE  3a.  Limits  of  o  and  p  for  the  zone -transport  system 

Range  of  a 


Value  ofk  L=0  L=oo 
e  e 


Enriching  —  P  +  1 

Section  . 


Range  of  P 


Value  ofk  L=0  L=<io 
s  s 


Sign  Sign 

of  B  of  /Ll 

s 


Stripping 

Section 


k  <  ■  — 

-f  -  w 


-  rnv 


kW  +  £(l-i 


b.  Limits  of  a  and  p  for  the  zone-void  system 


Range  of  a 


Value  of  k  L„  =  0 
e 


L  =  « 
e 


Enriching 

Section 


^4 


vnj/ 

1  +  w 

vnl/ 

1  +  w 


Range  of  p 


Value  of  k 


L,  =  0 


L  =  00 
e 


Sign  Sign 

of  B  of  u 

s 


Stripping 

Section 


kW  f(l-k) 
W 
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Table  3,  Table  2,  and  algebraic  manipulation,  where  required,  reveals  that 
the  limits  of  a  and  P  are  the  same  for  the  zone-void  and  the  zone-transport  sys¬ 
tem  when  L,  and  L  are  infinite.  When  L  =  0  the  value  of  p  in  the  zone-void 
system  is  k/  6  while  in  the  zone -transport  system  P  =  1.  In  the  stripping  section 
of  both  systems  when  L  =  0  the  actual  weight  of  material  in  the  section  is  £. 

The  difference  in  the  value  of  p  when  L  =  0  is  due  to  the  difference  in  the  method 

s 

of  waste  collection.  In  the  zone-transport  system  the  waste  is  collected  when  the 

last  zone  is  completely  molten.  The  solute  concentration  in  the  waste  then  is  equal 

to  feed  concentration  when  L  =  0  and  therefore  a  must  equal  unity.  In  the  zone- 

void  system  the  waste  is  not  collected  until  a  weight  of  material  £  -  w  of  the  last 

zone  has  frozen  then  the  remaining  liquid,  w,  is  collected  as  waste.  When  L  =0 

s 

the  whole  column  is  the  last  zone  and  the  last  zone  is  assumed  to  freeze  by  a  nor¬ 
mal  freeze.  It  is  easily  shown  that  under  these  conditions  p  must  equal  k/  6. 

When  Lg  =  0,  a  =an///l+u)  in  the  zone-void  system  while  in  the  zone -transport 
system  a  =  1.  This  difference  in  the  limit  of  a  is  due  to  the  different  mechanisms 
of  product  collection  and  the  change  which  these  mechanisms  make  on  the  mathe¬ 
matical  definition  of  the  solute  concentration  in  the  product.  Within  the  limits  of 
the  mathematical  models  both  limiting  values  of  q  are  correct. 

The  first  zone-void  system  to  demonstrate  Pfann's  hypothesis  was  built 
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at  this  laboratory  for  the  purification  of  Sil^.  Buford  and  Starks  later  continued 

this  effort  for  this  laboratory  on  a  contractual  basis.  In  an  effort  to  eliminate  some 

of  the  experimental  deterrents  to  ultrapurification  inherent  in  the  zone-void  system 
19 

Moates  et  al  designed  a  new  continuous  zone  refining  system.  A  complete  de- 

20 

scription  of  this  system  has  been  presented  elsewhere.  This  system  which  we 

refer  to  as  the  matter  transport  system  uses  the  matter  transport  phenomenon 
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described  by  Pfann  ’  as  the  driving  force  for  material  flow  in  either  the  en¬ 
riching  or  stripping  section  of  the  system.  Material  movement  in  the  other  section 
is  accomplished  by  means  of  voids. 

In  his  description  of  the  zone-void  system  Pfann  always  considered  =  P£- 
In  the  zone-void  system  if  the  change  in  density  which  occurs  in  going  from  the 
solid  to  the  liquid  state  (or  vice  versa)  is  taken  into  account,  material  can  be  made 
to  flow  even  if  the  voids  are  reduced  to  zero  volume  in  one  of  the  arms.  In  such 
an  arm,  material  flow  is  accomplished  by  collecting  an  amount  of  material  per  pass 
that  is  equal  to  the  difference  in  volume  between  a  zone  of  material  in  the  cold  and 
hot  condition.  Figures  7  and  8  show  two  systems  which  utilize  this  principle. 

In  order  to  simplify  the  description  of  the  operation  of  these  systems,  assume 
Pg  >P£  and  k  <  1.  Then  in  Figure  7  the  enriching  section  is  on  top,  the  heaters 
travel  downward  and  matter  transport  moves  material  from  the  mid-feed  upward 
to  the  product  exit.  At  the  start  of  each  cycle  the  uppermost  heater  is  in  position 
at  the  top  of  the  section.  It  melts  the  solid  starting  at  the  product  sidearm  down 
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Figure  7.  Matter  Transport  System  with  Feed  Section  Located  in  the  Middle 


em  with  Feed  Located  Above 
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to  a  plane  where  a  weight  of  solid  material  h  has  melted.  The  liquid  expands 
overflowing  into  the  product  sidearm  and  down  into  the  product  container.  As  the 
heater  travels  down  the  column,  solid  freezes  at  the  level  of  the  sidearm.  When 
the  liquid  zone  reaches  the  mid-feed,  it  joins  its  liquid  length  with  the  mid-feed 
and  then  reciprocates.  The  liquid  zone  that  was  brought  down  contracts  on  freezing, 
drawing  up  into  the  column  a  volume  of  feed  equal  to  the  volume  of  product  which 
has  been  collected  in  that  pass.  Repetition  of  this  process  (that  is,  the  passage  of 
molten  zones)  provides  both  purification  by  the  freezing  process  and  continuous 
material  movement  via  the  matter  transport  phenomenon. 

In  Figure  8,  the  apparatus  is  arranged  in  the  same  manner  as  a  conventional 
zone-void  refiner.  The  enriching  and  stripping  sections  are  both  below  the  feed. 

In  the  matter  transport  section,  the  enriching  section,  the  heaters  travel  upward. 
Material  movement  in  this  section  is  downward  and  product  is  collected  at  the  bottom. 
Material  movement  is  produced  by  the  same  steps  used  to  describe  Figure  7.  The 
only  difference  is  that  when  the  heaters  start  a  pass,  the  bottom  heater  melts  in  a 
weight  of  solid  material  h,  as  before,  however  the  expanding  material  instead  of 
overflowing  directly,  as  in  Figure  7,  pushes  the  molten  material  below  it  down  and 
then  up  over  the  inverted  U-tube  in  the  product  collector. 

The  stripping  sections  in  both  Figures  7  and  8  are  identical  and  use  examples 

g 

of  the  void-former  already  described. 

A  derivation  of  steady- state  equations  describing  the  enriching  section  of  a 
matter  transport  system  is  not  as  clear-cut  as  in  the  zone-void  or  zone-transport 
systems.  In  the  zone-void  and  zone-transport  systems  experimental  product  col¬ 
lection  provides  a  concise  mathematical  definition  of  the  solute  concentration  ex¬ 
isting  in  the  product  (Eqs.  (27)  and  (29)  respectively).  The  best  description  of  the 
solute  concentration  in  the  product  of  a  matter  transport  system  is  that  it  lies 

somewhere  between  that  predicted  by  the  zone-void  and  the  zone-transport  equa- 
22 

tions.  The  conservative  estimate,  which  the  authors  prefer,  is  given  by  the  zone- 
transport  steady-state  equations.  This  derivation  has  already  been  given  (Appen¬ 
dix  G). 

The  stripping  section  of  the  matter  transport  system  uses  a  void-former  and 
if  freezing  in  the  terminal  zone  is  by  normal  freeze  the  derivation  of  the  steady- 
state  equations  is  the  same  as  the  zone-void  system  (A-.  .’.idix  D).  If  the  void- 
former  described  in  Reference  9  is  used,  freezing  in  the  final  zone  of  each  pass  is 
not  accurately  described  by  the  normal  freeze  mechanism.  The  actual  freezing 
process  is  somewhere  between  a  normal  freeze  and  a  uniform  freeze.  If  the 
terminal  zone  freezes  by  a  uniform  freeze,  the  derivation  of  the  steady- state  equa¬ 
tions  would  be  that  given  in  Appendix  H. 

20 

The  expressions  used  to  define  material  flow  in  the  matter  transport  system 

5 

were  those  originally  used  by  Pfann  defined  on  a  weight  basis.  To  make  terminology 
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consistent  within  this  report  the  expressions  used  to  define  the  zone- void  system. 
Table  2.  can  be  used  to  describe  the  matter  transport  system  noting  that  in  the 
matter  transport  stripping  section  h  -  f  is  the  same  as  w  in  Table  2. 

For  the  three  continuous  zone  refining  systems  described  thus  far,  the  zone- 
void,  the  zone-transport,  and  the  matter  transport  systems,  only  two  sets  of  steady- 
state  equations  exist  for  each  section.  These  equations  are.  however,  relatively 
complex.  An  evaluation  of  the  effect  of  variation  of  some  of  the  pertinent  param¬ 
eters  by  the  iterative  substitution  of  various  values  of  the  parameters  in  question 
is  a  time  consuming,  arduous  task.  By  nondimensionalizing  to  obtain  the  minimum 
number  of  parameters  it  is  feasible  with  the  aid  of  a  high  speed  digital  computer, 
to  show,  in  graphical  form,  the  effect  of  the  variation  of  any  of  the  physical  param¬ 
eters  on  the  enriching  and  stripping  sections  of  each  system.  The  nondimension- 
aliz'i'^ion  transforms  the  steady- state  equations  of  each  section  into  a  pair  of 
equa.>ons  we  call  universal  equations.  An  equivalence  table  for  the  nondimensional 
variables  is  given  in  Table  4.  The  derivation  of  the  universal  equations  for  each 
section  is  given  in  Appendices  I-L. 


TABLE  4.  Nondimensional  terminology  equivalence 
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The  universal  equations  for  the  zone-void  enriching  section  are: 


vy 


k 


(37) 
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and 


k_^ 
1  -  y 


1  + 


1  -  k 


(vz  \ 


(38) 


For  the  stripping  section  of  the  zone-void  system  when  the  freezing  in  the  final 
zone  is  assumed  to  proceed  by  normal  freezing  the  universal  equations  are: 


k 


(39) 


and 


P  = 


e'^^d-ky) 


(l-ky)(l-y)^'*^  +  e'^^(l-y)  -  (1-y) 


(40) 


In  the  zone-transport  enriching  section  the  universal  equations  are  given  by 


-VJL 


e''  -  1 


(37) 


and 


a  = 


_ JL_LiS _ 

e'^^(l-k)  +  y  -  1 


(41) 


The  universal  equations  for  the  stripping  section  of  the  zone-transport  system 
or  the  zone-void  or  matter  transport  systems  when  the  terminal  zone  undergoes 
a  uniform  freeze  are 
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The  four  sets  of  universal  equations.  Eqs.  (37)  and  (38),  Eqs.  (39)  and  (40). 

£k]S.  (37)  and  (41)  and  Eqs.  (39)  and  (42)  were  programmed  and  run  on  a  high  speed 
digital  computer.  The  computations  were  made  for  k  >  1  as  well  as  for  k  <  1 . 

The  y  and  z  vadues  were  chosen  to  give  maximum  range  to  the  data.  Using 
Eqs.  (37)  and  (39)  the  families  of  curves  which  result  from  plotting  k  versus  v 
for  severad  y  vadues  au*e  shown  in  Figures  9  auid  10,  respectively.  For  plots  of 
a  or  p  versus  k  each  y  vadue  used  generated  a  family  of  z  curves.  There  were 
seven  y  vadues  used  for  each  of  the  four  cases.  For  k  <  1  the  twenty-eight  faunilies 
of  z  curves  which  resulted  are  shown  in  Figures  11-38  inclusive.  Figures  11-17 
show  a  versus  k  when  Eq.  (38)  is  used.  Figures  18-24  are  plots  of  P  versus  k 
when  Eq.  (40)  is  run.  Figures  25-31  au:e  the  q  versus  k  plots  of  Eq.  (41)  amd  Fig¬ 
ures  32-38  show  the  p  versus  k  curves  when  Eq.  (42)  is  used. 

A  compaurison  of  Figures  11-17  with  Figures  25-31  reveads  that  the  effect  of 
the  difference  in  the  definition  of  the  solute  concentration  in  the  product  on  the 
vadue  of  a  increases  with  decreasing  z.  However,  for  z  >  5  the  v^ue  of  a  for 
both  cases  is  almost  identical. 

A  comparison  of  Figures  18-24  with  Figures  32-38  permits  am  evaluation  of 
the  effect  of  assuming  a  normal  or  uniform  freeze  in  the  finad  zone  of  a  stripping 
section.  It  is  readily  seen  that  for  z  =  2  the  normal  freeze  section  provides  a 
slightly  better  sepauration  thaui  the  uniform  freeze,  however,  for  all  other  values 
of  z  the  values  of  P  aire  essentiadly  identical. 

The  functions  of  the  respective  sections  of  a  zone-refining  system  when  k  <  1 
are  reversed  when  k  >  1;  thus,  the  stripping  section  becomes  the  enriching  sec¬ 
tion  and  vice  versa.  Graphs  of  k  versus  a  or  P  for  each  set  of  universal  equa¬ 
tions  when  k  >  1  are  shown  in  Figures  39-66.  Thus: 

1.  Figures  39-45  are  plots  of  a  versus  k  when  Eq.  (38)  is  used. 

2.  Figures  46-52  are  plots  of  P  versus  k  when  Eq.  (40)  is  used. 

3.  Figures  53-59  are  plots  of  a  versus  k  when  Eq.  (41)  is  used. 

4.  Figures  60-66  are  plots  of  P  versus  k  when  Eq.  (42)  is  used. 

A  comparison  similar  to  those  made  with  Figures  11-38  can  be  made  with  Fig¬ 
ures  39-66. 
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Figure  51.  Plot  of  log  p  vs  k  when  Eq.  (40)  is  Used  and  k  >  1 
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Figure  57.  Plot  of  log  a  vs  k  when  Eq. 
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3.  RESULTS 


As  a  result  of  the  effort  described  in  this  report  several  accomplishments 
have  been  made  in  the  field  of  continuous  zone  refining.  An  error  was  found  in 
the  overall  solute  balance  equation  used  in  the  derivation  of  the  basic  equations 
describing  the  stripping  section  of  the  continuous  zone- void  system  reported  by 
Pfann.  The  correct  solute  balance  equation  was  reported  (Eq,  (20))  and  a  new 
set  of  steady-state  equations  derived  (Appendix  D).  A  difference  between  the  zone- 
void  and  the  zone- transport  system  not  previously  noted  was  found.  This  differ¬ 
ence,  which  is  in  the  mathematical  definition  of  the  solute  concentrations  of  the 
product  and  waste  of  each  system,  and  the  error  made  by  Pfami  in  the  solute  bal¬ 
ance  equation  of  the  stripping  section  of  the  zone- void  system  were  used  to  com¬ 
pletely  explain  the  discrepancies  between  the  limits  of  a  and  |3  obtained  by  using 
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Pfann's  steady-state  equations  for  the  zone-transport  system  and  Abe's  steady- 

state  equations  for  the  zone -transport  system.  An  error  in  Abe's  zone-transport 

equations  due  to  improperly  accounting  for  the  density  change  due  to  melting  was 

noted  and  corrected. 

A  new  continuous  zone  refining  system,  the  matter  transport  system,  was 
20 

described.  The  matter  transport  system  increases  the  flexibility  of  design  of 
a  zone  refining  system  which  can  result  in  increased  product  purity.  Steady-state 
equations  describing  the  three  continuous  systems  treated  herein  were  rederived 
with  particular  attention  being  paid  to  terminology  and  the  treatment  of  the  change 
in  density  due  to  melting.  The  feasibility  of  describing  the  matter  transport  sys¬ 
tem  by  the  mathematical  model  of  either  the  zone-void  or  the  zone -transport  sys¬ 
tem  was  noted,  however,  the  authors  have  indicated  their  preference  for  the  zone- 
transport  model  wherein  the  mode  of  freezing  in  the  terminal  zone  of  the  stripping 
section  is  taken  into  account. 

The  actual  values  of  a  and  P  obtained  using  the  matter  transport  system  will 
be  somewhere  between  those  predicted  by  the  zone-void  and  the  zone -transport 
systems. 

The  choice  of  the  optimum  physical  parameters  of  a  continuous  zone-refining 
system  required  to  produce  any  given  a  or  P  has  been  greatly  simplified  by  the 
introduction  of  nondimensional  parameters.  The  nondimens ionalization  permits 
complete  description  of  a  zone  refining  system  with  a  minimum  number  of  param¬ 
eters.  A  high  speed  digital  computer  can  and  has  been  used  to  obtain  results 
which  in  graphical  form  clearly  and  simply  indicate  the  magnitude  of  the  effect 
produced  by  variation  of  any  of  these  pareuneters.  The  nondimensionalization 
transforms  tho  steady- state  equations  of  each  section  of  a  continuous  system  into 
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a  pair  of  equations  which  we  have  called  universal  equations.  The  universal  equa¬ 
tions  for  the  zone-void  system  and  the  zone-transport  system  were  derived 
(Appendices  I-L).  The  universal  charts  which  result  from  computer  computations 
made  by  substituting  realistic  parameter  values  into  the  equations  are  shown  in 
Figures  9-66.  With  the  aid  of  these  charts  it  can  be  seen  that  the  difference  be¬ 
tween  the  enriching  section  of  the  zone-void  and  zone -transport  systems  is  negli¬ 
gible  for  y  <0.9  and  z  >  10.  In  the  stripping  section  the  difference  between  the 
normal  and  uniform  freeze  is  also  negligible  for  y  <  0.9  and  z  >  10,  however,  of 
more  importance  is  the  fact  that  the  curves  clearly  demonstrate  the  futility  of  in¬ 
creasing  z  over  10. 
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Appendix  A 

DERIVATION  OF  THE  SOLl'TE  DISTRIBLTION  IN  A  SEMI-INFINITE  INGOT 
AFTER  THE  PASSAGE  OF  ONE  MOLTEN  ZONE.  EQ.  (3)  (AFTER  PFANN^) 


The  purpose  if  this  derivation  is  to  find  an  expression  which  will  define  the 
solute  concentration  in  the  solid.  C,  at  any  plane  along  a  charge  after  one  molten 
zone  has  passed.  The  authors  have  modified  Pfa  ui's  original  derivation  to  account 
for  the  change  in  density  due  to  melting.  The  matter  transport  phenomenon  and  its 
complicating  effects  on  the  shape  of  the  charge  can  be  accounted  for  by  the  change 
due  to  melting. 

In  order  to  derive  an  equation  for  the  solute  distribution  under  these  conditions 
the  spacial  orientation  of  the  charge  and  the  ratio  of  the  density  of  the  solid  to  the  density 
of  the  liquid,  Pg/p£  .  must  be  considered.  If  Figure  1  is  placed  in  a  vertical 
position,  the  mathematical  treatment  is  fairly  straightforward.  If  k  <  1  and 
Ps  >  p£  a  molten  zone  can  be  made  to  travel  downwaud  without  fear  of  breakage 
due  to  matter  transport,  since  expansion  is  taken  care  of  in  the  first  zone.  When 
the  molten  zone  reaches  the  bottom  of  the  charge  the  zone  contracts  on  freezing. 
The  equation  for  the  solute  concentration  C  at  the  plane  which  exists  at  any  weight 
of  solid  material  x  measured  from  the  beginning  of  the  frozen  chi  .  ;  o  the  freez¬ 
ing  interface  of  the  molten  zone  being  passed  through  the  charge  is  derived  as 
follows: 

Let: 

=  the  average  solute  concentration  at  any  plane  in  the  charge, 
prior  to  passage  of  the  molten  zone,  expressed  in  weight 
fraction. 


A2 


Cjj  =  the  solute  concentration  in  the  molten  zone  in  weight- fraction, 
s  =  the  weight  of  solute  in  a  zone  at  any  x 
=  the  weight  of  solute  in  the  zone  at  x  =  0 
k  =  effective  distribution  coefficient 
£  =  weight  of  liquid  in  the  molten  zone 
Assume: 

1.  £  is  constant 

2.  the  charge  is  of  unit  cross-section 

When  the  molten  zone  moves  a  distance  equivalent  to  the  weight  increment  dx, 
a  weight  of  solid,  dx,  will  be  frozen  out  of  the  zone  at  the  freezing  interface  and 
an  equal  weight  of  solid  will  be  melted  in.  The  weight  of  solute  frozen  out  is 

kC^  dx  (A-1) 


but 


--  f  (A- 2) 

The  weight  of  solute  melting  into  the  zone  is 

C^dx  (A-3) 

Therefore  the  resulting  change  in  s  is 

ds  =  (c^-^)dx  {A-4) 


or 


ds  +  ^  dx  =  C  dx 
£  o 


Integrating  Eq.  (A-5)  yields 


kx 

£ 


s  “  T—  C  e 
k  o 


kx 

f 


(A-5) 


(A-6) 


Evaluating  Eq.  (A-6)  between  the  limits  0  and  x 


A3 


C 


. 

C  (k-1)  e  +  C 
o  o 


or 


(1-k) 


(A-7) 


(A-8) 


(A-9) 


(A-10) 


(3) 


Appendix  B 

DERIVATION  OF  THE  ULTIMATE  DISTRIBUTION  IN  A  BATCH  SYSTEM. 

EQS.  (4)  -  (6)  (AFTER  PFANn2) 


Consider  Figure  1;  assume  that  a  sufficient  number  of  passes  have  been  passed 
through  the  charge  to  establish  the  ultimate  distribution.  Under  these  conditions 
when  a  molten  zone  is  passed  through  the  ingot  it  does  not  change  the  solute  distri¬ 
bution.  Using  the  assumptions  and  definitions  of  terms  as  given  in  Appendix  A,  the 
solute  concentration  in  the  molten,  zone,  C^ix),  is  defined  as  C(x)/k.  C|(x)  is  also 
given  by  Eq.  (B-1).  Thus 

x+f 

Cj(x)=^  /  C(x)dx  (B-1) 


or 


C(x)  =  / 


x+£ 


C(x)  dx 


(B-2) 


The  solution  of  Eq.  (B-2)  is 
C(x)  =  Ae®* 


(4) 


substituting  Eq.  (4)  in  Eq.  (B-2)  and  integrating  yields 


B2 


AeBx.|keBX(e«f.l) 


(B-3) 


Solving  Eq.  (B-4)  for  x  --  0  gives 


e«^-l 


If  the  whole  ingot  is  considered  then  Eq.  (B-2)  can  be  written 


1 

C(x)  =  —  J  C(x)  dx 


(B-4) 


From  Eqs.  (4)  and  (B-4),  Eq.  (6)  is  readily  obtained  thus 


1 


Appendix  C 

DERIVATION  OF  THE  STEADY-STATE  EQI  ATIONS  FOR  THE  ENRICHING 
SECTION  OF  THE  ZONE-VOID  SYSTEM.  EQS.  (12)  -  (IS)  (AFTER  PFANN’) 


When  the  enriching  section  of  the  zone-void  system  has  reached  the  steady 
state,  removal  of  product  at  the  start  of  each  pass  form  a  void.  Passing  the  void 
through  the  section  shifts  the  solute  distribution  C(x)  as  shown  in  Figure  4.  When 
a  molten  zone  is  passed  through  the  section  it  shifts  the  solute  distribution  back 
to  the  steady  state. 

Under  these  conditions  the  derivation  proceeds  as  follows: 

Let; 

£  =  the  weight  of  liquid  material  in  the  molten  zone  that  is  maintained 
molten  and  caused  to  move  by  a  moving  heater. 

h  -£  =  weight  of  product  collected  per  pass. 

h  =  the  weight  of  material  equal  to  the  weight  of  liquid  material 

that  is  maintained  molten  and  caused  to  move  by  a  moving  heater 
plus  the  weight  of  product  collected  per  pass. 

Lg  =  weight  of  material  in  the  column  up  to  the  plane  where  the  solute 
concentration  is  kCj. 

The  amount  of  solute  in  a  molten  zone  having  its  freezing  interface  at  the 
plane  x  is  C(x)£/k.  This  amount  of  solute  is  equal  to  the  same  amount  of  solute 
that  was  melted  into  the  molten  zone  minus  thi  amount  that  had  frozen  out.  Thus 


C2 


The  solution  of  Eq.  (C’l)  is 


B  X 

C(x)  =  Aj  +  ® 


(C-2) 


Substituting  Eq.  (C-2)  into  Eq.  (C-1)  and  integrating  yields 


1  -e 


B  (h-1) 
e 


^2® 


B  X 
e 


/  B  h  \ 


(C-3) 


There  are  two  sets  of  terms  on  both  sides  of  Eq.  (C-3).  One  set  of  terms  ccpt?iins 

B  X 

the  variable  in  the  form  e  ;the  other  set  contains  constamt  terms  only,  in  order 

B  X 

for  the  equation  to  hold,  the  e  ®  terms  have  to  be  equal  and  the  constaint  terms 
must  be  equal.  _ 

£3  X 

•  »  0 

Equating  the  e  terms  and  simplifying  yields 


B  h  B  t 

e  =  1  + 

k 


(13) 


when  the  constant  terms  are  equated  and  rearranged 

0  =  A2W  - 


where 


(i)  = 


Bgf  (1-k) 


r  B^(h-f)i 
1  -  e 


(C-4) 


(14) 


B  L 

Using  Eq.  (C-4)  auid  the  assumed  boundary  condition  C(L  )  =  A,  +  A„e  ®  ® 

0  X 

expressions  for  and  A2  cam  be  determined.  Thus 


kCjOj 


B  L 
e  e  . 

e  +  oj 


kC 


f 


B  L 
e  e  . 

e  +  OJ 


(C-5) 


(C-6) 


C3 


The  expression  for  the  solute  concentration  in  the  product  is  given  by: 


h-l  /  B  x\ 

Cp(h-J)  -  f  '  jdx 


(C-7) 


Integration  of  Eq.  (C-7)  shows 


B  (h-l) 
e  ®  -  1 


(C-8) 


From  Eqs.  (C-5),  (C-6),  and  (C-8)  it  can  be  shown  that 


B  L 
,  e  e 


(C-9) 


or 


^e  "  0.434  B  ‘  ^ 

e  ' 


(12) 


where 


^  - 


kh  -  1 
h  -  £ 


(15) 


Appendix  D 

DERIVATION  OF  THE  CORRECTED  STEADY -STATE  EQUATIONS  FOR  THE 
STRIPPING  SECTION  OF  THE  ZONE-VOID,  SYSTEM.  EQS.  (21)  -  (24) 


The  steady-state  solute  distribution  C(x)  for  the  stripping  section  of  the  zone- 
void  system  is  shown  in  Figure  5.  The  shift  of  the  solute  distribution  due  to  pass¬ 
ing  a  void  through  the  system  is  indicated  by  the  dashed  curve.  The  travel  of  the 
molten  zone  through  the  system  shifts  the  solute  distribution  back  to  the  steady 

state.  C(x)  is  made  up  of  two  parts  joined  at  x  =  L  ,  the  zone  refining  portion 

s 

and  the  normal  freezing  portion. 

The  derivation  of  the  steady-state  equations  which  define  such  a  system  is  as 
follows;  Let 

w  =  weight  of  waste  collected  per  pass 

£  =  weight  of  liquid  material  in  the  molten  zone  that  is  maintained 
molten  and  caused  to  travel  by  a  moving  heater 

^s 

h  =  a  weight  of  material  equal  to  £  — 

P£ 

In  the  zone  refining  portion  of  C(x)  the  amount  of  solute  in  the  molten  zone, 
whose  freezing  interface  is  at  x,  is  given  by: 


D2 


^  --  «c  / 


x+l-w 


,  +  J  C(x)  dx  -  J  C(x)  dx 

£-w  0 


(D-2) 


The  solution  of  Eq.  (D-2)  is 


C  =  +  A^e 


(D-3) 


Substituting  Eq.  (D-3)  into  Eq.  (D-2)  and  integrating  yields 


V 


I  A  T  B  (f-w)]  A2  B  xj-  B  (£-w)  1  ( 

j'-T'  L'  -‘J 


(D-4) 


Equating  the  e  terms  of  Eq.  (D-4)  and  simplifying  gives 


B  (£-w)  B  £ 

s  ,  s 

e  =  — 


Equating  the  constant  terms  of  Eq.  (D-4)  and  using  Eq.  (22)  yields 


kCf  -  A^  +  A2 


(D-5) 


Equation  (D-5)  is  a  boundary  condition  describing  the  solute  concentration  at 
X  =  0.  In  order  to  solve  for  A^  and  A2,  a  second  expression  for  A^  and  A2  is 
needed.  The  solute  balance  over  the  entire  cycle  is  such  an  expression.  Thus; 


£C  =  /  C(x)  dx  +  wpc 

0  * 


(D-6) 


Using  Eq.  (D-3),  Eq.  (D-6)  can  be  rewritten 


.  £-Jw1  ^ 

f  £  -  w 


B  <£-w) 
s 


2  r 

i - 1  e  -  1 


(D-7) 


From  Eqs.  (D-7)  and  (D-5)  are  obtained 


A  -  kC, pP 
1  f 


(D-8) 


D3 


A2  =  kCjd  -  /3m) 


(D-9) 


where 


£  -  k(£-w) 


It  has  been  assumed  that  Eq.  (D-3)  describes  the  solute  concentration  in  the 
charge  up  to  the  plane  at  weight  of  material  L^.  The  solute  concentration  in  the 
molten  zone  from  which  the  solid  at  is  freezing  has  the  solute  concentration 
C(L  )/k.  Assuming  that  between  x  =  L  and  x  =  L  +  £  a  normal  freeze  takes 
place,  then  the  normal  freezing  equation  (£q.  (2))  describes  the  solute  concentration 
over  this  distance.  Figure  5.  Since  a  weight  of  material  equal  to  w  is  removed 
from  the  waste  end  of  the  section  in  order  to  form  a  void  during  each  pass,  then 
the  solute  concentration  in  the  plane  +  £  -  w  is  equal  to  The  solute 

concentration  in  the  solid  at  the  plane  +  £  -  w  is  kCj.  Therefore  since: 


C(Lg)  =  Aj  + 


(D-10) 


and  the  normal  freezing  equation  is 


C  =  kC^(l-g)" 


combining  Eqs.  (D-10)  and  (2)  gives 


^w  =  PCf  =  -k^  V' 


(D-11) 


C(Lg)  =  p6Cj 


(D-12) 


where 


‘  •  '■(f) 


D4 


Substituting  Eqs.  (D-8).  (D-9).  and  (D-12)  in  Eq.  (D-10)  gives 


e 


®s^s  _  p(6-ku) 
■  itU-PM) 


or 


L 

s 


1 

0.434  B 

s 


klTTO'J 


(D-13) 


(21) 


I  . 

I 

i 
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Appendix  E 

DERIVATION  OF  THE  STEADY-STATE  EQl  ATIONS  FOR  THE  ENRICHINO 
SECTION  OF  THE  ZONE-TRANSPORT  SYSTEM  tSING  ABE'S*^  TECHMQI  E 


Abe  points  out  that  when  the  steady  state  has  been  obtained,  if  we  consider 
the  two  zones  a  and  b  (Figure El),  the  total  amount  of  solute  in  these  two  zones 
when  a  is  liquid  and  b  is  solid  is  equal  to  the  total  amount  of  solute  in  the  zones 
when  a  is  solid  and  b  is  liquid.  This  equality  assumes  that  the  shape  of  the  ingot 
is  a  rectangular  parallelepiped  and  that  surface  tension  is  negligible.  This  expres¬ 
sion  is  written 

+  /  C(x+P+£)dx=  /  ^  C(x)  dx  +  (E-1) 

X  X 

or 


.x+P+£ 


C(x+P+£)  -  C(x)  =  j  J  [C(x+P+£)  -  C(x)  dx] 


(E-2) 


The  rate  of  change  of  solute  in  the  molten  zone  is: 


d  =  [C(x+P+£)  =  C(x)]  dx 


(E-3) 


E3 

The  solution  of  Eq.  (E-2>  is 

B  X 

C(x+P+f  >  -  C(x>  =  Ae  ®  (E-4) 

Substituting  Eq.  (E-4>  into  Eq.  (E-2)  and  simplifying  yields 
B  (P+£>  B  £ 

e  ®  =  k  1  (33) 

Equation  (E-3)  can  be  rewritten  thus 

^  [C(x+P+Jf)  -  C(x)J  (E-5> 


Substituting  Eq.  (E-4)  into  Eq.  (E-S)  and  integrating  between  the  limits  0  and  x 
yields 


C(x)  =  C(o) + 


(E-6) 


Assuming  the  boundary  condition C(Lg)=  kC^  and  using  Eq.  (E-6)  it  can  be  shown 
that 


/  B  £ 

‘-e  ■  I;  I*  *  -  '^<‘>>1  i 

By  definition 

1  r'’*' 

‘Cp-  ph-f  C(x)  dx 

From  Eqs.  (E-8)  and  (E-6)  it  can  be  shown  that 

S  ^  ^  B  (P+£) 

e  e 


(E-7) 


(E-8) 


(E-9) 


Consideration  of  the  gain  and  loss  of  material  in  the  two  regions  0  <  x  <  P  +  I, 
and  P+£  <  x  <  2(P+£)  when  they  are  molten  and  solid  respectively  yields 


E4 


^  2(P+f) 

Cl  +  /  C(x)dx  = 

P  p*f 


/  C(x)dx  + 
0 


k 


(E-10) 


Using  Eq.  (E-8)  and  performing  the  necessary  integration  and  simplification  it  can 
be  shown  that  Eq.  (E-10)  can  be  rewritten  as 


r  +  A  -  C(P^f) 

p  k  ■  k 


(E-11) 


When  X  -  0,  Eq.  (E-4)  becomes 

C(P+£)  -  C(o)  =  A  (E-12) 

Substituting  Eq.  (E-12)  into  Eq,  (E-11),  using  Eq.  (E-9)  to  substitute  for  C(o)  and 
solving  for  A  yields 


A  = 


aCjB^I(i-k)(P+£) 
i  -  k<P+£) 


Using  Eq.  (E-9),  Eq.  (E-7)  can  be  rewritten 


^  ■  r  *" 

e 


B  £ 


(kC,-C  )  + 


Ak  '""'f'-p'  '  k(P+f) 


From  Eq.  (E-13) 


t  -  k(P+£) 
kaCj(l-k)(P+£) 


(E-13) 


(E-14> 


(E-15) 


From  Eqs.  (E-15)  and  (E-14)  it  can  be  shown  that 
j  /[f-k(P+£)]  (k-Q) 


e 


k(P+£)\  £a(l-k) 


+  1 


or 


B  L 
.  e  e 


j  /[£  -k(P^£)]  (k-g) 


k(P+£)  \ 


fa(l-k) 


+  1 


(E-16) 


(E-17) 


Appendix  F 

OEftlVATION  OF  THE  STEADY-STATE  EQUATIONS  FOR  THE  STRIPPING 
SECTION  OF  THE  ZONE-TR  ANSPORT  SYSTEM  I  SING  ARE’S'2  TECHNIQUE 


As  in  Appendix  E,  consider  the  amount  of  solute  in  the  two  zones  c  and  d,  of 
Figure  FI,  when  c  is  liquid  and  d  is  solid  and  when  c  is  solid  amd  d  is  liquid. 
Assume  the  ingot  shape  is  a  parallelopiped  except  in  the  molten  zone  where  the 
shape  of  the  liquid  zone  is  as  shown.  Then  the  solute  balance  can  be  expressed  as 


C(x)(h+w) 

— c — 


x+h 


x+h 


+  f  C(x+h)  dx  =  J  C(x)  dx  + 


(h+w)  C{x+h) 


(F-l) 


or 


C(x+h)  -C(x)  = 


/  (c(x+h)-C(x)^  dx 


(F-2) 


For  the  molten  zone  the  rate  of  change  of  solute  is 


^^C(x)[h^\  ^  /c(x+h)- C(x)\  dx 


(F-3) 


The  solution  of  Eq.  (F-2)  is 


F3 


C(x+h)-C(x)  =  Ae 


B  X 
s 


From  Eqs.  (F-4)  and  (F-2)  it  can  be  shown  that 


B^h  B  (h+w) 
e  ®  =  - -  +  1 


Equation  (F-3)  can  be  rewritten  thus 


^  {C(x+h)-c<x)) 


Substituting  Eq.  (F-4)  in  (F-6)  and  integrating  between  0  and  x  yields 


C(x)  =  C(o)  + 


Ak 


B  (h+w) 
s 


e  -  1 


At  X  =  0  the  following  boundary  condition  holds 


C(o)  =  kC, 


By  definition 


Cw  = 


C(^) 


and 


s  s 


From  Eqs.  (F-7),  (F-8)  and  (F-9)  it  can  be  shown  that 
C, 


Cw  =  ^  =  PC,  .  C,  * 


E  L 

e  ®  ®-l 


Solving  Eq.  (F-11)  for  L 


(F-4) 


(F-5) 


(F-6) 


(F-7) 


(F-8) 


(F-9) 


(F-10) 


(F-11) 


F4 


L 


s 


/  B  Cf(B-l)(h+w) 


(F-12) 


For  the  two  zones  L  =  1  <  x  <  and  L  <  x  <  L'  the  gain  and  loss  of  mate- 
s  s  s  s  “ 

rial  when  they  are  molten  and  solid  respectively  can  be  expressed  as 


-Ls  (h+w)C(j^ 

wpCf  +  J  C(x)  dx  =  - ; - ■*’ - 

^  L  -h 
s 


(F-13) 


Using  Eqs.  (F-5)  and  (F-7)  and  performing  the  indicated  integration  Eq.  (F-13) 
becomes 


I  (‘ -  h^)  ■  S 


solving  for  A  yields 


B^Cj(h+w)(w[  1-p]  +h[  1-k] ) 
w  +  h(l-k) 


Substituting  Eq.  (F-15)  in  Eq.  (F-12)  and  simplifying  yields 


-  ±  r  Ph(l-k) 
s  ■  h(l-k)  +  w(l-p) 


Ph(l-k) 

h(l-k)  +  w(l-p) 


(F-14) 


(F-15) 


(F-16) 


e 


(F-17) 


Appendix  G 

DERIVATION  OF  THE  STEADY-STATE  EQLATIONS  FOR  THE  ENRIEIIING 
SECTION  OF  THE  ZONE-TRANSPORT  SYSTEM  USING  PFANN’S'^  TECHNIQIIE 


This  derivation  uses  the  same  mathematical  technique  shown  in  Appendix  C. 
The  modifications  are  imposed  by  the  variation  in  boundary  conditions  and  termi¬ 
nology.  Thus  the  amount  of  solute  in  the  molten  zone  whose  freezing  interface  is 
at  the  plsuie  x  is  C(xK/k  auid  is  given  by  the  solute  balance  expression,  thus 


=  /  C(x)  dx  -  /  C(x)  dx  -  5^  /  C(x)  dx 

o  o  o 


(G-1) 


The  solution  of  Eq.  (G-1)  is 


C  =  Aj  +  A2  e 


(G-2) 


Substituting  Eq.  (G-2)  in  Eq.  (G-1)  and  integrating  yields 


f  w  v\  vr^2 


n  ^  B^:<  ^ 

.  Be  ^ 


A„  BjP+i) 
— ^  e  ® 

Be 


A„f  A„f  B^(P+£) 

B  (P+£)  B  (P+£)  ® 
e  e 


(G-3) 


G2 


Equating  the  e  terms  Eq.  (G-3)  gives 


B  (P+f)  B  !. 


Equating  the  constant  terms  in  Eq.  (G-3)  yields 


(G-4) 


where 


One  of  the  boundary  conditions  is 


kCj  =  Aj  +  A2  e 


(G-S) 


From  Eqs.  (G-4)  and(G-5)it  can  be  shown  that 


kC,w 
'  “b~l 

e  e  , 

e  +  a> 


(G-6) 


'2  B  L 
e  e 

e  +  u> 


(G-7) 


The  solute  concentration  in  the  product  is 


,  P+£  /  B  x\ 

PT7  /  (^A^  +  A2e  « 


(G-8) 


Integrating  Eq.  (G-8)  and  utilizing  Eqs.  (G-6),  (G-7)  and  (33)  to  solve  for  yields 


®e^e  _  u)k(P+£)  +  £ 

®  "  q(P+£) - 


(G-9) 


Appendix  H 

DERIVATION  OF  THE  STEADY-STATE  EQUATIONS  FOR  THE  STRIPPING 
SECTION  OF  THE  ZONE-TRANSPORT  SYSTEM  USING  PFANN’S''*  TECHNIQUE 


In  the  zone -transport  system  the  waste  is  withdrawn  before  the  terminal  zone 
begins  to  freeze.  Since  back  reflection  is  not  considered  then,  when  the  density 
change  due  to  melting  is  taken  into  accoimt,  the  derivation,  which  is  very  similar 
to  that  given  in  Appendix  D,  is  as  follows: 

The  amount  of  solute  in  a  molten  zone  having  its  freezing  interface  at  x 
(Figure  HI)  is 


fcTx)  /■* 

— r —  =  fCj  +  j  C(x-[n-£+w])  dx-  j  C(x)  dx 

h  o 


Equation  (H-1)  can  be  rewritten  as 


(H-1) 


C(x)  + 


x+Cw 


fCj  +  J  C(x)dx-  J  C(x)  dx 
f -w  o 


(H-2) 


The  solution  of  Eq.  (H-2)  is 


B  X 


=  +  Age 


C 


(H-3) 


Weight  of  Solid  Material  (x) 


Figure  HI.  Steady-State  Solute  Distribution  in  the  Stripp:.  Section  of  the 
Zone  Transport  System  when  the  Terminal  Zone  Undergo.)^  a  Uniform  Freeze 


H3 


Substituting  Eq.  (H-3)  in  Eq.  (H-2)  and  integrating  gives 


B„x  , 

A  +  A  e  ®  =  — 

^1  ^2®  £ 

^2 

')1 

s 

V  = 

/J 

(H-4) 


B  X 
s 


Equating  the  e  terms  yields 


B  (f-w)  B  f 
e  ®  1+  -S- 

k 


(H-5) 


Equating  the  constant  terms  of  Eq.  (H-4)  and  using  Eq.  (H-b)  gives 


kCf  -  Aj^  +  A2 


(H-6) 


the  boundary  condition  existing  at  x  =  o.  The  solute  balance  for  the  entire  cycle  is 


.£-w 


£Cj  =  J  C(x)dx  +  wpC 


f 


(H-7) 


From  Eqs.  (H-3)  and  (H-7)  it  cam  be  shown  that 

B„(£-w)  ^ 


'f 


£  -  Pw 


A,  + 


A,  (e 


1 "  b;t£^i 


(H-8) 


Using  Eqs.  (H-6)  amd  (H-8)  it  is  readily  shown  that 


A^  =  ’  Cjp)i 


(H-9) 


and 


A2  =  kCj(l-pfi) 


(H-10) 


where 


w 


^  "  i.  -  k(£-w) 


(H-11) 


H3 


Substituting  Eq.  (H-3)  in  Eq.  (H-2)  and  integrating  gives 


B  X  ,  Ao  /  B 

Ai  +  Aje  "  =1 


A„  /  B„(£-w)  A,  B^x/  B^(f-w) 


"2  “s* / 

«  V 


(H-4) 


Equating  the  e  terms  yields 


B  (1-w) 

e  ®  1  +  -f 


(H-5) 


Equating  the  constant  terms  of  Eq.  (H-4)  and  using  Eq.  (H-b)  gives 


kCj  =  +  Ag 


(H-6) 


the  boundary  condition  existing  at  x  =  o.  The  solute  balance  for  the  entire  cycle  is 


£Cj  =  /  C(x)dx  +  wpC^ 


(H-7) 


From  Eqs.  (H-3)  and  (H-7)  it  can  be  shown  that 


P  £-pw  .  .  .  ^2  ^ 

^f  f -  w  ■  B  (i-w) 

s 


B  (£-w) 
s 


(H-8) 


Using  Eqs.  (H-6)  and  (H-8)  it  is  readily  shown  that 


A^  =  ’  CjpM 


(H-9) 


A2  =  kCj(l-pfi) 


(H-10) 


where 


(H-11) 


-13)  it  is  straightforward  to  show  that 


Appendix  I 

DERIVATION  OE  THE  UNIVERSAL  EQUATIONS  FOR  THE  ENRICHING 
SECTION  OF  THE  ZONE-VOID  SYSTEM  (EQS.  (37)  and  (38)) 


The  pertinent  equations  are 


®e^e  /li/  .V 
e 


B  h  B  £ 

e  =  1  +  -r- 
k 


Bg£(l-k) 


1  -  e 


B^(h-£) 


4,  = 


kh  -  £ 
h  -  £ 


where 


a 


(C-9) 

(13) 

(14) 

(15) 


rewriting  Eqs.  (C-9),  (13),  (14)  and  (15) 


12 


e®*'  =  1  +  Bh  • 

h  k 


“(-i) 


Bh.  r 
n 


1  L 

Using  the  equalities  Bh  =  v,  ^  =  y,  =  z  to  rewrite  Eqs.  (I-D-  (1-4)  yields 


respectively 


■  ■) 


e'  .  .  .  f 


Substituting  Eqs.  (1-7)  and  (1-8)  in  Eq.  (1-9)  and  solving  for  a  gives  Eq.  (38)  thus 


rearranging  Eq.  (1-6)  Eq.  (37). 


e'^  -  1 


Appendix  J 

DERIVATION  OF  THE  UNIVERSAL  EQUATIONS  FOR  THE  STRIPPING 

SECTION  OF  THE  ZONE-VOID  SYSTEM 


The  basic  equations  when  the  terminal  zone  freezes 


®s^s  .  |l(6-ku) 
■  k(l-PM) 


B  (£-w)  B  £ 

e«  =1.-^ 


by  a  normal  freeze  are 

(D-12) 

(22) 

(23) 


and 


^  "  £  -  k(£-w) 


(24) 


where 


J2 


Letting  Bg£  =  v,  y  =  .  and  z  =  -p  Eqs.  (D-12),  (22).  (23)  and  (24)  can  be 

rewritten  respectively  as 


6  =  kd-y)^"*' 


SuLstituting  Eqs.  (J-4)  and  (J-3)  in  Eq.  (J-1)  and  solving  for  p  gives 


e^^(l-ky) 


(l-ky)(l-y)^'‘'  +  e''^(l-y)-(l-y) 


rearranging  Eq.  (J-2)  yields 


(J-1) 

(J-2) 

(J-3) 

(J-4) 

(40) 

(39) 


Appendix  K 

DERIVATIO^  OF  THE  LISIVERSAL  EQUATIONS  FOR  THE  ENRICHING 
SECTION  OF  THE  ZONE-TRANSPORT  SYSTEM 


The  pertinent  equations  are 

.®e^e  .  a)k(P+£)  +  £ 

®  -  a  iP^t) - 

B  (P+f)  B  £ 

e  ®  =  1  +  — f— 

k 

and 

,  P 

(P+£)(l-k) 

£ 

Letting  Bg(P+£)  =  v,  p^rj  ~  V  pTf  '  (G-9).  (33)  and  (32)  can  be 

rewritten  respectively 

^vz  ^  ^ 

a 

f  +  1  (K-2) 


(G-9) 

(33) 

(32) 


K2 


(K-3) 


Substituting  Eq.  (K-3)  into  Eq.  (K-1)  and  solving  for  a  yields 


y  -  k 


e'^^(l-k)  +  y  -  1 


rewriting  Eq.  (K-2)  gives 


e  -  1 


Deriving  Eqs.  (41)  and  (37)  from  the  basic  equations  which  result  when  the  zone- 
transport  steady-state  equations  are  derived  by  Abe's  technique  (Appendix  E)  can 
be  accomplished  as  follows. 


B  L 
.  e  e 


k(P+£) 


j  f  [£-k(Pt£)]  (k-g) 


£a(l-k) 


(E-17) 


B  (P+£) 


=  1  + 


V  — 0 

When,  as  above,  B^(P+£)  -  v,  "  y  pTl  ’  (E-17) can  be  rewritten 


•{(yfei") 


(K-4) 


Solving  Eq.  (K-4)  for  o  yields  Eq.  (41).  Thus 


y  -  k 

e'^^(l-k)  +  y  -  1 


It  was  shown  above  that  Eq.  (33)  can  be  transformed  into  Eq.  (37).  Therefore, 
when  carried  to  the  universal  equation  form,  the  identity  of  the  results  of  Abe's 
and  Pfann's  derivation  techniques  is  obvious. 


I 


I 


■9---*rrs 


Appendix  L 

DERIVATION  OF  THE  UNIVERSAL  EQUATIONS  FOR  THE  STRIPPING 
SECTION  OF  THE  ZONE-TRANSPORT  SYSTEM  (EQS.  (39)  AND  (42)) 


The  basic  equations  are 

_®s^s  .  p(l-fi) 

1-pfi 

(H-15) 

B  (£-w)  B^£ 

e  -  1  +  — 

(35) 

^  ■  kCTrrr 

(36) 

Letting  =  v. 


vz  _  P(1-m) 
■ 


vy  1  V 
e  =  1  +  r 
k 


.  1  -  y 

1  -  ky 


c 

=-j^  =  y  and  -|p  =  z,  Eqs.  (11-15),  (35)  and  (36)  can  be  rewritten 

(L-1) 

(L-2) 

(L-3) 


I  1 

[ 

I 


Substituting  Eq.  (L-3)  in  Eq.  (L-1)  and  solving  for  p  yields 
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